Test of Lorentz invariance with spin precession of ultracold neutrons 
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A clock comparison experiment, analyzing the ratio of spin precession frequencies of stored ul- 
tracold neutrons and 199 Hg atoms is reported. No daily variation of this ratio could be found, 
from which is set an upper limit on the Lorentz invariance violating cosmic anisotropy field 
b± < 2 x 10~ 20 eV (95% C.L.). This is the first limit for the free neutron. This result is also 
interpreted as a direct limit on the gravitational dipole moment of the neutron \g„\ < 0.3 eV/c 2 m 
from a spin-dependent interaction with the Sun. Analyzing the gravitational interaction with the 
Earth, based on previous data, yields a more stringent limit |g n | < 3 x 10 -4 eV/c 2 m. 

PACS numbers: 14.20.Dh, 11.30.Er, 11.30.Cp, 28.20.-v 

Lorentz symmetry is a fundamental hypothesis of our 
current understanding of physics and is central to the 
foundations of the Standard Model of particle physics 
(SM). However, the SM is widely believed to be only 
the low energy limit of some more fundamental theory, 
a theory which will probably violate more symmetries 
than the SM, in order to accomodate some features of 
the universe currently lacking in the SM, e.g., the baryon 
asymmetry. A SM extension including Lorentz and CPT 
violating terms has been presented in [if. It provides a 
parametrisation of effects suitable to be tested by low 
energy precision experiments. In particular, clock com- 
parison experiments 0, S| have proven to be particularly 
sensitive to spin-dependent effects arising from a so-called 
cosmic spin anisotropy field b filling the whole universe. 
This Letter reports on a search for such an exotic field 
via its coupling to free neutrons. 

In the presence of a field b, the two spin states of the 
neutron will encounter an extra contribution to the en- 
ergy splitting corresponding to the potential V = a ■ b 
where a are the Pauli matrices. Thus, if a neutron is sub- 



jected to both a static magnetic field B and the new field 
b, its spin will precess at the modified Larmor frequency 
/„, which to first order in b is given by 

fn = ^ B + h-*. (i) 

Jn 2tt h B w 

We searched for a sidereal modulation (at a period of 
23.934 hours) of the neutron Larmor frequency induced 
by b± , the component of b orthogonal to the Earth's ro- 
tation axis. The experiment is also sensitive to a possi- 
ble influence of the Sun on the spin precession dynamics, 
leading to a solar modulation (at a period of 24 h) of 
the Larmor frequency, as proposed in [4{. Such an effect 
could arise from a non-standard spin-dependent compo- 
nent of gravity |5L M or from another long-range spin- 
dependent force [7], [8|. In particular, a non-zero neutron 
gravitational dipole moment g n would induce a coupling 
through (see also |9() 

Vgdm = 511—3- cr • r, (2) 
where G is Newton's constant, and for the mass M and 
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the distance r we use the Sun mass M Q and the distance 
Earth-Sun tq . 

The experimental apparatus at the PF2 [l(| beam- 
line at ILL, Grenoble, is normally used to search for the 
electric dipole moment of the neutron (nEDM) [ll], [13] • 
The apparatus permits spin-polarized ultracold neutrons 
(UCN) to be filled into a volume, stored, and then 
counted and classified by polarization state. While con- 
fined, the neutrons can be exposed to static (normally 
B w 1 f/T), as well as to oscillating, magnetic fields. A 
surrounding four layer Mu-metal shield suppresses the 
external magnetic field and its fluctuations. Although b 
acts like a magnetic field influencing the particles' spin 
precession, it can, per definition, not be suppressed by 
the Mu-metal shield. 

The neutron Larmor frequency / n — j n B/(2ir) w 
30 Hz, is measured via the Ramsey method of separated 
oscillatory fields 13, ljj]. Following filling, an initial oscil- 
lating field pulse rotates the neutron spin by 7r/2, leaving 
the magnetic moment at right angles to the static holding 
field B, whereupon it precesses. Following a free spin pre- 
cession time of typically 100 s, a second oscillating field 
pulse, phase coherent with the first pulse, further rotates 
the neutron spin by it/ 2. The accumulated phase is mea- 
sured by counting the populations of the two resulting 
spin states following the second Ramsey pulse. For each 
cycle about 10 4 neutrons are counted allowing a measure- 
ment of / n with a statistical precision of A/ n w 50 /iHz. 
A unique feature of this nEDM apparatus is the use of 

Within the neutron 
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a mercury co-magnetometer 

precession chamber, nuclear spin polarized 199 Hg atoms 
precess in the same B field as the neutrons. The Larmor 
frequency /n g = 7H g -B/(27r) ss 8 Hz is measured opti- 
cally for each cycle to a precision of A/n g ~ 1 /iHz. The 
pumping and analysing light are generated by two lamps 
filled with 204 Hg and Ar plasma. In addition, four scalar 
Cs magnetometers [HI are placed above and below the 
precession chamber (see Fig. [T]). They provide on-line 
measurements of the magnetic field with a precision of 
150 fT and were used to measure the vertical gradients 
of the magnetic field. 

For the clock comparison experiment, we use the ra- 
tio R = f n /fHg which suppresses the effect of magnetic 
field fluctuations in the limit of a perfectly homogeneous 
field. The existing constraints for 199 Hg [2] are suffi- 
ciently tight, so within the sensitivity of this experiment, 
new physics effects can only show up in f n . While the 
Earth is rotating, together with the vertical quantization 
axis, the new physics effects under consideration would 
appear in a harmonic change of R: 
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(3) 



The constant term 5R would be induced by the com- 
ponent of a new field parallel to the Earth rotation axis, 
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FIG. 1: Vertical cut through the cylindrical precession cham- 
ber for UCN and 199 Hg. Schematically indicated are the 
~ 70 mm diameter Cs vapor filled bulbs and their mounts. 
The scalar Cs magnetometer measures the magnitude of B 
found at the center of the spherical bulb. 



while the amplitude A would be induced by the orthog- 
onal component. Since magnetic field gradients related 
effects can be mistaken for the steady term SR, we first 
focus on the search for a non zero amplitude A. 

Data were recorded in April-May 2008 with the B field 
pointing downwards. The first 35 h of data were recorded 
starting on April 21 07h20 UT, followed by a break of 
255 h, and then 85 h of uninterrupted data were col- 
lected. As we search for a signal modulation in R(t), the 
runs were combined after subtracting the mean values R 
of the corresponding runs. Figure [2] shows folded data, 
and its discrete spectral analysis. The error bars indi- 
cate combined statistical errors of the neutron and the 
Hg frequency extraction, dominated by the former one. 
The spectral analysis shows that no particular frequency 
can be extracted from the data and the whole dataset is 
compatible with a signal of null amplitude (Xnuii = 0.98). 

The neutron frequency extraction fl4| requires a fit 
of the visibility a of the Ramsey resonance curve. The 
value of this parameter depends on the value of the mag- 
netic field gradients. In order to avoid systematic ef- 
fects correlated with the value of these gradients (which 
could be daily modulated), a was fitted in small (typi- 
cally one hour) subsets of data to ensure that the neutron 
frequency extraction does not create any bias. 

The R parameter also depends on the value of the mag- 
netic field gradients inside the chamber: while the center 
of mass of the thermal 199 Hg gas coincides with the cham- 
ber center, the UCN center of mass is about h = 2.8 mm 
lower, due to gravity. This offset is related to the vertical 
gradient dB/dz [121 ]: 
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Daily variations of dB/dz would be the main system- 
atic uncertainty in our analysis and could appear, e.g., 
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FIG. 2: The upper figure shows the variation of the ratio 
R around its average. For clarity, the data are folded mod- 
ulo 24 hours and binned every half hour. The lower figure 
shows modulus of the discrete Fourier Transform of the same 
dataset. the line would be an evidence that this frequency is 
too much represented as compared with a white noise. 



due to a daily modulation of the Earth magnetic field. 
Therefore, the vertical gradients were monitored by the 
Cs magnetometers (two on top of the storage chamber 
and two below it). At the frequency of interest, 1/24 h, 
the amplitude of fluctuations of dB/dz was mesured to 
be < 20pT/m resulting in a daily modulation of R with 
an amplitude < 2 x 10~ 7 , according to Eq. Q. This 
effect is small enough to prevent the R ratio departing 
from a white noise signal, as one can see in Fig. [2] Other 
possible sources of false daily modulated signal have been 
investigated and ruled out. Besides magnetic field inho- 
mogeneities, the main remaining effect is related to the 
light shift of the mercury frequency. We estimated the 
associated relative error to be « 10 -7 with our analyzing 
light intensity. The drifts in intensity of the light was 
measured to be less than 10 %. 

To extract an upper limit for the daily modulation am- 
plitude, a frequentist confidence level analysis was per- 
formed on the unfolded data. The method consists in 
determining wether a given signal hypothesis (a given 
amplitude A and phase 4>) can be excluded at 95% C.L. 
when compared with the null hypothesis. This method 
is known to optimally discriminate two signal hypotheses 
[l6j |. For a given A and (/>, we form the quantity: 
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(5) 



R[i] - Asm(2irt\i]/T - 
AR[i\ 



where N = 2070 is the total number of data points, R[i] 
is the R ratio subtracted from its mean value in indi- 
vidual datasets and T is either the solar period or the 
sidereal period. From the measured data, Qdata(A 4>) is 
calculated. We consider the probability distribution of Q, 
Pnuii(Q) in the null hypothesis and p s i gna ,i(Q) in the sig- 
nal hypothesis. These two probability distributions have 
been calculated (for each signal hypothesis) using Monte- 
Carlo simulations. The confidence level of the hypothesis 
(A, cf>) is then defined as 
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The amplitude A is excluded at 95% C.L. if, for all phases 
4>, we have CL(A,(f>) < 0.05. The statistical limit ob- 
tained this way is A < 0.58 x 10~ 6 , both, for the sidereal 
and solar period. Accounting for the systematics that 
dB I dz modulations could counteract any other modula- 
tion, the amplitude due to new physics is limited to: 



A < 0.8 x 10" 6 95% C.L. 



(7) 



Our result can be interpreted in terms of a limit on 
the cosmic spin anisotropy field for the free neutron. In 
this case, T is the sidereal period and the amplitude A is 
related to the b± component according to: 



A = 2b 



cos(A) 



(8) 



where A = 45°12'22" is the latitude of the experiment in 
Grenoble and h is the Planck constant; thus: 



b , < 2 x 10" 20 eV 95% C.L. 



(9) 



Table |T] compares this result to existing limits on other 
particles. The result reported here is the first limit for 
the free neutron. It is complementary to the more precise 
atomic experiments 0, S| that can be interpreted as lim- 
its concerning bound neutrons inside nuclei. Contrary to 
the results of@, Q , the neutron result is free from model- 
dependent nuclear corrections and possible related sup- 
pression effects. Being a factor 200 more stringent than 
the limit for the proton, it is the best free nucleon limit 
to date. 

The result Eq. ([7]) can also be interpreted as a limit 
on the gravitational dipole moment of the neutron. In 
this case, T is the solar period and the amplitude A is 
expressed as 



A = 2 9n 



GMq cos(A) 



(10) 



where the inclination of the Earth with respect to the 
ecliptic plane, supressing the effect by less than 5%, is 
neglected. We thus obtain the following upper bound 



\g n \ < 0.3 eV/c 2 m 95% C.L. 



(11) 



4 



Reference 


System 


Particle 


b ± [eV] 


Berglund et ai, [2] 


Hg 


& Cs 


bound neutron 


9 x 10" 










electron 


2 x 10" 


20 


Bear et al., [3] 


Xe 


& He 


bound neutron 


2 x 10" 


22 


Phillips et al..\17] 




H 


proton 


4 x 10" 


18 


Heckel et ai, [181 




e 


electron 


7 x 10" 


22 


Bennet et al, [19] 






positive muon 


2 x 10" 


15 








negative muon 


3 x 10" 


15 


This analysis 


n 


fcHg 


free neutron 


2 x 10" 


20 



TABLE I: Results of more restricting upper limits (at 
95% C.L.) on b±(e), b±(N), b±(p), b±(n), b±(n) the couplings 
between a cosmic spin anisotropy field and different particles. 



In principle, much more stringent limits can be set using 
the Earth as a source of the new spin-dependent effect. 
In this case, one has to search for a steady signal contri- 
bution 5R added to or substracted from the main cou- 
pling term 7 n /7H g , see Eq. ([3]), depending on the direc- 
tion of the magnetic field. Previous measurements with 
the same apparatus constrained the difference of the val- 
ues Ro, the value of R for dB/dz = 0, for magnetic field 
pointing upwards and downwards to T3, 23| 



\SR\ = - \Rqi - Roi \ < 1.6 x 1(T 6 



95% C.L. (12) 



Using the Earth's mass (M s ) and radius (r^) in Eq. ([2]), 
this can be translated into a limit on the neutron gravi- 
tational dipole moment 



\9n\ = \SR\ 



2 GMa 



< 2.5 x 1(T 4 eV/c 2 m 95% C.L. 

(13) 

Finally, the limits derived above can be discussed in 
terms of the Hari-Dass framework of spin dependent 
gravity [f|: 

^Hari-Dass = »iGM + a * GM A ( 14 ) 
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where a\ and ot2 are dimensionless parameters and v is 
the neutron velocity with respect to the source. ol\ is 
directly proportional to the neutron gravitationnal dipole 
moment and the limit (1131) leads to: 
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2\g n \c 2 
he 



< 2.5 x 10 3 95% C.L. 



(15) 



this is the best limit for the neutron although still far 
above the natural value a± « 1. A more stringent limit 
|q;i| <2x 10 2 was obtained using 199 Hg and 201 Hg [21(, 
however involving nuclear model uncertainties. 

While ai is best constrained using the Earth as a 
source, the search for daily modulation is the natural 



way to probe 6*2, using the Sun as the source. Our limit 
(fTTj) translates to 

|a 2 |<3xl0 10 95% C.L. (16) 
As for the limit on the cosmic anisotropy field, this is the 
first limit on the free neutron. 

We are grateful to the ILL staff for providing us 
with excellent running conditions and in particular ac- 
knowledge the outstanding support of T. Brenner. We 
also benefitted from the technical support through- 
out the collaboration. We acknowledge fruitful dis- 
cussions with R. Lehnert. This work is supported by 
the French-Polish funding agencies (French: convention 
IN2P3-Laboratoires Polonais, grant No. 05-120 and Pol- 
ish: COPIN/283/2006), the original development of the 
nEDM apparatus was funded by the UK's PPARC (now 
STFC) and the Cs magnetometers were financed by the 
Swiss National Science Foundation (200020-119820). 



* On leave of absence from INRNE, Sofia, Bulgaria 
t On leave from PNPI, Russia 

* Also at University of Zurich 
§ pignol@lpsc.in2p3.fr 

■ rebreyend@lpsc.in2p3.fr 
** roccia@lpsc.in2p3.fr 

[1] D. Colladay and V. A. Kostelecky, Phys. Rev. D 55, 6760 
(1997). 

C. J. Berglund et al, Phys. Rev. Lett. 75, 1879 (1995). 

D. Bear et al, Phys. Rev. Lett. 85, 5038 (2000). 
R. Golub, Inst. Phys. Conf. Ser. No. 42, 104 (1978). 
J. Leitner and S. Okubo, Phys. Rev. 136, B1542 (1964). 
N. D. Hari Dass, Phys. Rev. Lett. 36, 393 (1976). 

B. A. Dobrescu and I. Mocioiu, JHEP 0611, 005 (2006). 
J. E. Moody and F. Wilczek, Phys. Rev. D 30, 130 
(1984). 

R. Golub, D. J. Richardson, S. K. Lamoreaux, Ultra-cold 
Neutrons, (Adam Hilger, Bistol, 1991). 

A. Steyerl et al, Phys. Lett. A 116, 347 (1986). 
P. G. Harris et al, Nucl. Inst. Method, in Phys. Res. A 
440, 479 (2000). 

C. A. Baker et al, Phys. Rev. Lett. 97, 131801 (2006). 
N. F. Ramsey, Phys. Rev. 76, 996 (1949); 78, 695 (1950). 
K. Green et al, Nucl. Inst. Meth. in Phys. Res. A 404, 
381 (1997). 

S. Groeger, A. S. Pazgalevl, A. Weis, Appl. Phys B 80, 
6 (2005). 

C. Amsler et al. [Particle Data Group], Phys. Lett. B 
667, 1 (2008). 

D. F. Phillips et al, Phys. Rev. D 63, 111101(R) (2001). 

B. R. Heckel et ai, Phys. Rev. Lett. 97, 021603 (2006). 
G. W. Bennett et al. [Muon (g-2) Collaboration], Phys. 
Rev. Lett. 100, 091602 (2008). 

C. A. Baker et ai, Phys. Rev. Lett. 98, 149102 (2007). 
B. J. Venema et al, Phys. Rev. Lett. 68, 135 (1992). 



[2] 
[3] 
K 
[5 
[6] 
[7] 



[9 

[10 
[11 

[12 
[13 
[14 

[is; 

[16 

[17 
[18 
[19 

[20 
[21 



